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Abstract 

A function J with simple and nice algebraic properties is defined on a subset 
of the space of complex sequences. Some special functions are expressible in 
terms of first of all the Bessel functions of first kind. A compact formula in 
terms of the function J is given for the determinant of a Jacobi matrix. Further 
we focus on the particular class of Jacobi matrices of odd dimension whose 
parallels to the diagonal are constant and whose diagonal depends linearly on 
the index. A formula is derived for the characteristic function. Yet another 
formula is presented in which the characteristic function is expressed in terms 
of the function J in a simple and compact manner. A special basis is constructed 
in which the Jacobi matrix becomes a sum of a diagonal matrix and a rank-one 
matrix operator. A vector-valued function on the complex plain is constructed 
having the property that its values on spectral points of the Jacobi matrix are 
equal to corresponding eigenvectors. 
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1 Introduction 

The results of the current paper are related to the eigenvalue problem for finite- 
dimensional symmetric tridiagonal (Jacobi) matrices. Notably, the eigenvalue problem 
for finite Jacobi matrices is solvable explicitly in terms of generalized hypergeometric 
series p]. Here we focus on a very particular class of Jacobi matrices which makes 



it possible to derive some expressions in a comparatively simple and compact form. 
We do not aim at all, however, at a complete solution of the eigenvalue problem. 
We restrict ourselves to derivation of several explicit formulas, first of all that for 
the characteristic function, as explained in more detail below. We also develop some 
auxiliary notions which may be, to our opinion, of independent interest. 

First, we introduce a function, called 5, defined on a subset of the space of complex 
sequences. In the remainder of the paper it is intensively used in various formulas. 
The function $ has remarkably simple and nice algebraic properties. Among others, 
with the aid of $ one can relate an infinite continued fraction to any sequence from 
the definition domain on which $ takes a nonzero value. This may be compared to the 
fact that there exists a correspondence between infinite Jacobi matrices and infinite 
continued fractions, as explained in [2] Chp. 1]. Let us also note that some special 
functions are expressible in terms of First of all this concerns the Bessel functions 
of first kind. We examine the relationship between J and the Bessel functions and 
provide some supplementary details on it. 

Further we introduce an infinite antisymmetric matrix, with entries indexed by 
integers, such that its every row or column obeys a second-order difference equation 
which is very well known from the theory of Bessel functions. With the aid of function 
5 one derives a general formula for entries of this matrix. The matrix also plays an 
essential role in the remainder of the paper. 

As an application we present a comparatively simple formula for the determinant of 
a Jacobi matrix of odd dimension under the assumption that the neighboring parallels 
to the diagonal are constant. As far as the determinant is concerned this condition 
is not very restrictive since a Jacobi matrix can be written as a product of another 
Jacobi matrix with all units on the neighboring parallels which is sandwiched with 
two diagonal matrices. The formula further simplifies in the particular case when the 
diagonal is antisymmetric (with respect to its center). In that case zero is always an 
eigenvalue and we give an explicit formula for the corresponding eigenvector. 

Finally we focus on the rather particular class of Jacobi matrices of odd dimension 
whose parallels to the diagonal are constant and whose diagonal depends linearly 
on the index. Within this class it suffices to consider matrices whose diagonal is, 
in addition, antisymmetric. In this case we derive a formula for the characteristic 
function. Yet another formula is presented in which the characteristic function is 
expressed in terms of the function $ in a very simple and compact manner. Moreover, 
we construct a basis in which the Jacobi matrix becomes a sum of a diagonal matrix 
and a rank-one matrix operator. This form is rather suitable for various computations. 
Particularly, one can readily derive a formula for the resolvent. In addition, a vector- 
valued function on the complex plain is constructed having the property that its values 
on spectral points of the Jacobi matrix are equal to corresponding eigenvectors. 

2 The function $ 

We introduce a function $ defined on a subset of the linear space formed by all complex 
sequences x = {xk}^ =1 - 
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Definition 1. Define $ : D — > 



oo oo 



#0) = 1 + X^ _1 ) m Yl Yl ■■■ Y x kl x kl+1 x k2 x k2+1 ...x km x km+1 (1) 

m=l fci=l fc 2 =fci+2 fc m =fc m _i+2 

where 

{oo 
{x k }^ =1 ; ^ |^fe^fc+i| < oo 
k=l 

For a finite number of complex variables we identify $(xx, X2, ■ ■ ■ , x n ) with $(x) where 
x = (xi, X2, ■ ■ ■ , x n , 0, 0, 0, ... ). By convention, we also put 5(0) = 1 where is the 
empty sequence. 

Remark 2. Note that the domain D is not a linear space. One has, however, £ 2 (N) C D. 
To see that the series on the RHS of ([I]) converges absolutely whenever x G D observe 
that the absolute value of the mth summand is majorized by the expression 



CO 

XjXj + i | 



^ I I 1 (f 

/ y \X kl X kl+ iX k2 X k2+ i . . . x km x km+ i\ < — I 2 J 

k&\ m \j=l 
kx<k 2 <...<k 

I'll 

Obviously, if all but finitely many elements of a sequence x are zeroes then $(x) 
reduces to a finite sum. Thus 

= 1, $(xi, x 2 ) = 1 - X X X 2 , $(Xi, x 2 , x 3 ) = 1 - x x x 2 - x 2 x 3 , 
$(xi, x 2 , x 3 , x A ) — 1 - x x x 2 - x 2 x 3 - x 3 x A + xiX 2 x 3 x A , etc. 

Let T denote the truncation operator from the left defined on the space of all 
sequences: 

T({x k }T =1 ) = {xk + i}T =l . 

T n , n = 0, 1, 2, . . . , stands for a power of T. Hence T n ({x k } k % 1 ) = {^fc+n}^!- 
The proof of the following proposition is immediate. 

Proposition 3. For all x G D one has 

$(x) = $(Tx)-x 1 x 2 d(T 2 x). (2) 

Particularly, if n > 2 then 

${Xi,X 2 , X 3 ,..., X n ) = $(x 2 , X 3 ,..., X n ) - XiX 2 $(x 3 , . . . ,X n ). (3) 

Remark 4. Clearly, given that #(0) = $(xi) = 1, relation ([3]) determines recursively 
and unambiguously $(xi, . . . , x n ) for any finite number of variables n G Z + (including 
n = 0). 
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Remark 5. One readily verifies that 

$(xi,x 2 , ...,x n )= $(x n , . . .,x 2 ,x 1 ). (4) 
Hence equality (jHJ) implies, again for n > 2, 

^(^l, . . . , X n —2, Xn—li x n) = St^-l; • • • j ^n-2; ^n-l) %n-l%n ^{Xl, • • • , £71-2)- (5) 



00 
k=0 



Remark 6. For a given x E D such that #(a;) 7^ let us introduce sequences {Pk} 
and {Qfc}£L by P = and P fe = #(:T2, . . . , for > 1, Q k = ■ ■ ■ ,x k ) for 

k > 0. According to ([5]), the both sequences obey the difference equation 

Y k+1 = Y k - x k x k+ iY k -i, k = 1, 2, 3, ... , 

with the initial conditions Po = 0, Pi = 1, Qo = Qi = 1, and define the infinite 
continued fraction 

^ - lim * - 



2^3 



1 _ X 3 X 4 

1 - . . . 

Proposition [3] admits a generalization. 

Proposition 7. For ever?/ x £ D and fceN one /ias 

#0) = . . . , x k ) $(T h x) - . . . , Xk-i)x k x k +i $(T k+1 x). (6) 

Proof. Let us proceed by induction in k. For fc = 1, equality ([6]) coincides with ([2]). 
Suppose (J6j) is true for fceN. Applying Proposition [3] to the sequence T k x and using 
([5j) one finds that the RHS of ([6]) equals 

$(xi, ...,x k ) $(T k+1 x) - $(x!, x k )x k+1 x k+2 $(T k+2 x) 

- . . . , x k -!)x k x k+1 $(T k+1 x) 

= $(x lt ...,x k , x k+ i) $(T k+1 x) - x k )x k+1 x k+2 $(T k+2 x). 

This concludes the verification. □ 

Remark 8. With the aid of Proposition [3] one can rewrite equality <Q as follows 

$(x) = $(Xi, . . . , X k ) $ f ^fi ' ' ' ,Xk ~^ Xk Xk+ i, X k+2 , X k+3 , . . . J . (7) 

V . . . ,x k ) J 

Later on, we shall also need the following identity 
Lemma 9. For any n G N one has 

ui$(u 2 , u 3 ,..., U„)#(t>i, f 2 , t>3, • • • , u n ) - u 2 , u 3 ,..., u n )$(v 2 , v 3 ,..., v n ) 

= ^ I J^Mfcffc j (Uj ~ Vj) $(u j+1 ,Uj +2 , ■ ■ .,U n )$(Vj +1 ,V j+2 , ■ ■ .,V n ). (8) 

j=i Vfe=i / 
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Proof. The equality can be readily proved by induction in n with the aid of (EJ). □ 
Example 10. For t, w £ C, \t\ < 1, a simple computation leads to the equality 

v im(2m-l)„..2m 



m=l 



This function can be identified with a basic hypergeometric series (also called q- 
hypergeometric series) defined by 

ai; 0)k ■ ■ ■ (°n q) k (, ^k„ ik(k-i)\ 1+s ~ r zk 



_ (6i; g) fc ■ ■ ■ (b s ; q) k \ J (q;q) t 

where r,s£ Z + (nonnegative integers) and 



fc-i 
i=o 



see [5]. In fact, the RHS in (jUJ) equals o0i(; 0; t 2 , —tw 2 ) where 

°° qk(k-i) 00 qHk-i) 

o*fc <* * *> = E 21 = £ (i _ 9)(1 _ ^ . . . (i _ 

with g, z G C, |o| < 1, and the recursive rule ([2]) takes the form 

o<M; 0; 0, z) = <M; 0; o, gz) + z <M; °; a, <7 2 ^)- (10) 

Put e(g; z) = o0i(i 0; q, (1 — 0)2). Then lim^x e(g; z) = exp(z). Hence e(g; z) can be 
regarded as a q-deformed exponential function though this is not the standard choice 
(compare with [5] or [6] and references therein). Equality (flOj) can be interpreted as 
the discrete derivative 

e(g;z) - e(q;qz) 2 

t, v = e{q;q z). 

(1 - q)z 

Moreover, in view of Remark [HI one has 

1 = o0i(;O;g,^) 

o0i(;O;g,z) 



1 + ^ 



1+ q2z 



1 + ... 

This equality is related to the Roger s-Ramanuj an identities, see the discussion in [3J 
Chp. 7]. 
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Example 11. The Bessel functions of the first kind can be expressed in terms of 
function More precisely, for v ^ — N, one has 

The recurrence relation ([2]) transforms to the well known identity 

zJ v (z) - 20 + l)J„+i(z) + zJ u+2 {z) = 0. 
To prove (1111) one can proceed by induction in j = 0, 1, . . . , m — 1, to show that 

oo oo oo 

E E E 

fel=l fc 2 =A:i+2 fc m =fc m _i+2 

1 

(U + fci)(l/ + fci + 1)0 + fc 2 )0 + ^2 + 1) • • • + fcm)0 + &m + 1) 
j oo oo oo 

fci— 1 &2 — ^i+2 k rfl —j=k 7n —j — i-\-2, 

l 

X , 
+ fa)(v + fci + 1)0 + ^2)0 + A; 2 + 1) • • • + ^m-jOO + k m-j + 1) 

1 

v 

+ fc m _j + 2)0 + fc m _j + 3) . . . + &m-j + j + 1) ' 

In particular, for j = m — 1, the RHS equals 

1 00 1 
(m - 1)! ^ (1/ + fc x )0 + fa + 1)0 + fa + 2) . . . + fa + m) 



1 ro + i) 



ml (v + 1)0 + 2) . . . + 77i) ^i! T(z/ + m + 1) 

and so 



W I \ _ w 2m +1 , 



ro + i) Vl^ + ^J fc =i/ m!r(i/ + m + l)' 

as claimed. Furthermore, Remark [6] provides us with the infinite fraction 
v + 1 ^+1(2™) 1 



w J u (2w) w 



2 



+ 1)0 + 2) 



w 2 



+ 2)0 + 3) 

2 



x _ + 3)0 + 4) 



1 - 



This can be rewritten as 



Mz) 2(„ + l) 

2(u + 2) — 



2(i/ + 3) - 



Z 2 



2(i/ + 4)-... 



Comparing to Example [Hi one can also find the value of 5 on the truncated 
sequence {w/(v + fc)}J? =1 . 

Proposition 12. For n G Z + and i/gC \ {— n, —n + 1, . . . , —1} one has 

^ • • • = e (-!)• n W-*). 

(12) 



v i/ + l'i/ + 2'"''z/ + 7i/ r(z/ + n + l)^ v y s!(n-2s)! 
In particular, for m, n G Z +; m <n, one has 



\ , [("-"i)/2] / x, / x, 

w w «A = m! (n-s)l(n-m-s)! ^ 2s 

m + l'm + 2' 'ny n! s! (m + s)! (n — m — 2s)! 

Proof. Firstly, the equality 

^ (n + 1 - k){n + 2 - k) . . . {n + s - 1 - k) 
^ (v + k)(v + k + l) . . .(v + k + s) 

n (n + 1) . . . (n + s — 1) 
s {v + n + s) {y + \){y + 2) . . . (u + s) 



k=l 



(14) 



holds for all n G Z + , u E C, u ^ — N, and s G N. To show ( I14p one can proceed by 
induction in s. The case s = 1 is easy to verify. For the induction step from s — 1 to 
s, with s > 1, let us denote the LHS of (JT4"]) by Y s (i/,n). One observes that 

,., , i/ + n + s — 1 T , . . z/ I rx I 2 s 1 . „ . 

r s (z/, n) = Y s -.i{v, n) r s _i(i/ + 1, n). 

s s 

Applying the induction hypothesis the equality readily follows. 
Next one shows that 

n-2s+2 n-2s+A n 

E E ••• E 

fel=l ^2=^1+2 fc s =fc a _i+2 

1 , x 

x 

{y + fci)(i/ + fcj + l){y + k 2 ){y + k 2 + 1) . . . {y + k 8 ){y + k s + 1) 

(n - 2s + 2){n - 2s + 3) . . . (n - s + 1) 

s! (u + l)(v + 2) . . . (u + s) (u + n - s + 2)(u + n - s + 3) . . . (u + n + 1) 
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holds for all n G Z + , s G N, 2s < n + 2. To this end, we again proceed by induction in 
s. The case s — 1 is easy to verify. In the induction step from s — 1 to s, with s > 1, 
one applies the induction hypothesis to the LHS of ffT5]) and arrives at the expression 



n-2s+2 j 



E 

k=l 



(u + k)(v + k + l)(8-l)\ 



(n - k - 2s + 3)(n - k - 2s + 4) . . . (n - A; - s + 1) 



(i/ + + 2)(z/ + A; + 3) . . . (u + k + s) (u + n - s + 3)(u + n - s + 4) . . . (u + n + 1) ' 



Using (|T4"1) one obtains the RHS of (fT5|) . as claimed. 

Finally, to conclude the proof, it suffices to notice that 



s [n/2] n -2s+l n-2s+3 n-1 

*(^t.^.-.^) = ^E(-d- E E - E 

' s=l fci=l k 2 =ki+2 fc s =fc s _i+2 



29 



(i/ + fci)(i/ + fci + + fc 2 )(f + h + 1) . ■ • (v + /c s )0 + k a + 1) 
and to use equality fTToT) . □ 

One can complete Proposition [T2] with another relation to Bessel functions. 
Proposition 13. Form,n G Z +; m <n, one has 

Ti\ ( 11) 11) II) \ 

m! \m + 1 m + 2 ray 

- V (m-^-l)!(n-m + 2 a + l)! ^ + o(u , m+n+1 j H) 

^ s\{n + s + l)\{n-m + s + v v y/ 

Proof. Recall the following two facts from the theory of Bessel functions (see, for 
instance, [U Chapter VII]). Firstly, for /i, v — N, one has 

IZo s! + s + l)r(z/ + s + 1) V2/ 
where (a) s = a(a + 1) . . . (a + s — 1) is the Pochhammer symbol. Secondly, for n G Z+, 

7rY n (z) = |- ( J v (z) - (-l) n J_,(z)) 
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For m, n G Z + , m < n, a straightforward computation based on these facts yields 

\(n— m— l)/2l , . , , . . 

, (n- s - 1)! (n-m- s - 1)\ / Z \ m ~ n + 2 

TX.J rr) [z)Y n (z , , , 

: ..._ + s ! n-m-2s-l ! V2 

s=0 



1 J s!(m + s)!(n-m-2s-l)! V: 

Elm - s - 1)! (n - m + 2s)! /z\«~ m + 2s n _ , N _ . . . /z\ 
,!(„ + .)i(„- m + ,)l (2) + 2M*)U*)lo g (-) (17) 



s=0 

00 



^ ' s\(m + s)\!n + s)\!m + n + s)\\2J \ K ' 

s=0 

—ip{m + s + 1) — ^>(n + s + 1) — ^(m. + n + s + 1) — ip{s + 1) 

where ip(z) = T'(z)/T(z) is the digamma function. The proposition follows from (|17p 
and (USD. □ 



Remark 14. Note that the first term on the RHS of ( 116]) contains only negative powers 
of u>. One can extend ( ITS]) to the case n = m — 1. Then 



nJ m (2w)Y m (2w) = - V (m . 8 ^ ( ^ )! + 0( W 2m logH). 

' (s! r (m + sj! 



=0 



3 The matrix 3 

In this section we introduce an infinite matrix 3 that is basically determined by two 
simple properties - it is antisymmetric and its every row satisfies a second-order dif- 
ference equation known from the theory of Bessel functions. Of course, in that case 
every column of the matrix satisfies the difference equation as well. 

Lemma 15. Suppose w G C \ {0}. The dimension of the vector space formed by 
infinite-dimensional matrices A = {A(m, n)} m , n& i satisfying, for all m,nGZ, 

wA(m, n — 1) — nA(m, n) + wA(m, n + 1) = (18) 

and 

A(n,m) = —A(jn,n), (19) 

equals 1. Every such a matrix is unambiguously determined by the value A(0, 1), and 
one has 

\/n G Z, A(n, n + 1) = A(0, 1). (20) 

Proof. Suppose A solves ( 118]) and (fT9]) . Then A(m, m) = 0. Equating m = n in (TIB"]) 
and using f[T9"]) one finds that A(n, n + 1) = — A(n, n — 1) = A(n — 1, n). Hence (|2"U|) 
is fulfilled. Clearly, the matrix A is unambiguously determined by the second-order 
difference equation (ITS]) in rz and by the initial conditions A(m, m) = 0, A(m, m+1) = 
A(0, 1), when m. runs through Z. 
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Conversely, choose A G C, A 7^ 0. Let A be the unique matrix determined by (I18p 
and the initial conditions A(m,m) = 0, A(m,m + 1) = A. It suffices to show that A 
satisfies f|T9|) as well. Note that A(m,m — 1) = —A. Furthermore, 

wA(m — 1, m + 1) — mA(m, m + 1) + wA(m + 1, m + 1) 
= wA(m — 1, m + 1) — mA(m — 1, m) + wA(m — 1, m — 1) 
= 0. 

From (|T8j) and the initial conditions it follows that A(m,m + 2) = (m + i)A/u>, and 
so mA(m, m + 2) = (m + l)A(m — 1, m + 1). Consequently, 

«M(m — 1, m + 2) — mA(m, m + 2) + Wv4(m + 1, m + 2) 

= -uM(m — 1, m + 2) — {m + l)A(m — 1, m + 1) + a)A(m — 1, m) 

= 0. 

One observes that, for a given m G Z, the sequence 

x n = — y4(m — 1, n) H A(m, n), n£Z, 

solves the difference equation 

to„_i - nx n + wx n+ i = (21) 

with the initial conditions x m +i = A{m + 1, m + 1), a; m+ 2 = v4(m + 1, m + 2). By the 
uniqueness, x n = A(m + 1, n). This means that, for all m, n G Z, 

«M(m — 1, n) — mA(m, n) + u>A(m + 1, n) = 0. 

Put B(m,n) = —A(n,m). Then 13 fulfills ffT8|) and B(m,m) = 0, B(m,m + 1) = A. 
Whence 5 = A. □ 

Lemma 16. Suppose w G C \ {0}. // a matrix A = {A(m, n)} mi „ e z satisfies (TTSj) 
(ED i/ien 

Vm,n G Z, A(m, -n) = (-l)M(m,n), A(-m,n) = (-l) m A(m,n). (22) 

Proof. For any sequence {x n }n€Z satisfying the difference equation fT2TT) one can verify, 
by mathematical induction, that x~ n = (— l) n x n , n = 0, 1, 2, . . .. □ 

Definition 17. For a given parameter wgC \ {0} let 3 = {3(m, n)} mngZ denote the 
unique matrix satisfying (TT8]) . ( TT9]) and 5(m, m + 1) = 1, Vm G Z. 

Remark 18. Here are several particular entries of the matrix 

m + 1 (m + l)(m + 2) 

J(m, m) = 0, ^(m, m+1) = 1, .,j(m,m+2) = , ^(m,m+3) = - 1, 

w w z 

with m G Z. Some other particular values follow from ( II 9p and ( I22p . Below, in 
Proposition [221 we derive a general formula for 3(m, ^)- 
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Lemma 19. For < m < n one has (with the convention 

( n ~ X ) ! .„m-n+lj WW W 



3(m, n) = ± w m ~ n+1 y , , . . . , . (23) 

m! \m + l m + 2 n — 1/ 

Proof. The RHS of (1231) equals 1 for n = m + 1, and (m + l)/w for n = m + 2. 
Moreover, in view of (jHJ), the RHS satisfies the difference equation (I2ip in the index 
n. □ 

Remark 20. From (P23J and (HQ) it follows that 

VmeZ, lim — 7 3(m,n) = J m (2w). 

n-voo (n — 1)1 

This is in agreement with the well known fact that, for any w G C, the sequence 
{J n (2w)} ne z fulfills the second-order difference equation (|2"T|) . 

Remark 21. Rephrasing Proposition [T31 and Remark [141 one has, for m, n G Z + , m < n, 

j (n w , /n x N (ra — s — l)!(n — m + 2s)! „_ m+2s 

vrJ m 2w r n 2w = -w \j(m,n)- > v — -rr^w" m+2s 

s\ (n + s)\ (n — m + s)\ 

+ 0(w m+n logH). 

Since, by definition, the matrix 3 is antisymmetric it suffices to determine the 
values 3(m, n) for m < n, m,n G Z. In the derivations to follow as well as in the 
remainder of the paper we use the Newton symbol in the usual sense, i.e. for any 
z G C and n G Z + we put 

z\ z(z — 1) . . . (z — n + 1) 



Proposition 22. Form,n G Z, m < n ; one /ias 

[(n-m-l)/2] , 1 \ / iM 

3(m,n)= V (-!)'( V"~ m ~ S ~ 1)! W — (24) 

Proof. We distinguish several cases. First, consider the case < m < n. Then (|24|) 
follows from (|2"3"|) and (TT5|) . Observe also that for m = n, m, n G Z, the RHS of 
is an empty sum and so the both sides in (1241) are equal to 0. 

Second, consider the case m < < n. Put m = —k, k G Z + . The RHS of (J21 
becomes 

[(n+*-l)/2] / N , , 

y (-i)-f V s - 1 y^+fe-^-i)' w - fc -, +2 , +1 (25) 



N n + — 2s — 1/ s! 

Suppose k < n. Then the summands in (1251) vanish for s = 0, 1, . . . , k — 1, and so the 
sum equals 

'"' ' s+k (w-fc-s-1)! (n-s-1)! fc _„ +2s+1 



(n- A; -2s- l)!s! (s + k)\ 
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■ w 



By the first step, this expression is equal to (— 1) $(k,n) = $(—k,n) (see Lemma [To]) . 
Further, suppose k > n. Then the summands in fl25l) vanish for s = 0, 1, . . .,n — 1, 
and so the sum equals 

" " " J ' -s-1 \ (k-s-])\ k . 



V (-iW ~ ) 

^ v ; \k-n-2s-lj (n + s)\ 



w 



Using once more the first step, this expression is readily seen to be equal to 
(-l) k+1 3(n,k)=Z(-k,n). 

Finally, consider the case m < n < 0. Put m = —k, n = —£, k,£ G Z+. Hence 
< £ < k. The RHS of (El becomes 



a '"'" I,/21 J -t-s-1 \(k-£-s-l)\ ,_ k+2s+1 



E (-13 



s=0 



fc-^-2s-l7 s! 



Using again the first step, this expression is readily seen to be equal to (— l) k+e+1 ^(£, k) 
= 3{-k,-£). □ 

4 The characteristic function for the antisymmet- 
ric diagonal 

For a given d G Z + let E± denote the (2d + 1) x (2d + 1) matrix with units on the 
upper (lower) parallel to the diagonal and with all other entries equal to zero. Hence 

{E+) j>k = 8j+i,k, ( E -)j,k = fy,fc+i» h k = -d, -d + 1, -d + 2, . . . , d. 

For y = (y_ d , y_ d +i,V-d+2, ■ ■ ■ , Vd) G C 2d+1 let diag(y) denote the diagonal (2d + 1) x 
(2d + 1) matrix with the sequence y on the diagonal. Everywhere in what follows, / 
stands for a unit matrix. 

First a formula is presented for the determinant of a Jacobi matrix with a general 
diagonal but with constant neighboring parallels to the diagonal. As explained in the 
subsequent remark, however, this formula can be extended to the general case with 
the aid of a simple decomposition of the Jacobi matrix in question. 

Proposition 23. For d G N, w G C and y = (y_ d , y- d +i, V-d+2, ■ ■ ■ , Vd) £ C 2d+1 , 
Ylk^iVkV-k 0, one has 

det ( diag(y) + wE + + wE_) = FT y k y_ k y #( — — ) #( — ...,_) 

Vfc=i / L ^ m Vd ' y - d ' 

w 2 / w w \ w w \ w 2 / w w \ / w w 



l/i V2/2 i/d/ V2/-i y-d/ y-i Vyi l/d/ \y-2 



(26) 



12 



Proof. Let us proceed by induction in d. The case d = 1 is easy to verify. Put 
J\fd(w; y) = det ( di&g(y) + wE + +wE_) . Suppose ( 126]) is true for some c? > 1. For given 
io G C and ?/ G C 2d+3 consider the quantity A/rf+i(w; y). Let us split the corresponding 
[2d + 3) x (2d + 3) Jacobi matrix into four blocks by splitting the set of indices into 
two disjoint sets {—d — 1, d + 1} and {— d, —d + 1, — d + 2, . . . , d}. Applying the rule 



det = det(A) det(D - CA~ l B) 

one derives the recurrence relation 

J\f d+1 (w;y- d -i,y-d,y-d+i, ■ ■ -,yd,yd+i) = yd+iy-d-M(w;yi d ,y-d+i,y-d+2, 

where 



,yd-i,y d ) 



Vd 



l 



w 



Vd, y- d 



l - 



w 



y-d- 



VdVd+ij \ y-dV-d-i. 

Now it is sufficient to use the induction hypothesis jointly with the equality 



'1 - x n _ix n ) £ x%,x 2 



. Xn-2, 



, %2) ■ ■ ■ j 3^)1—1; -En) 



1 %n—l%n, 

which is valid for n > 2 and which follows from relations (j4j) and ([7]), with k = 2. □ 
Remark 24. Let us consider a general finite symmetric Jacobi matrix J of the form 

\ 



/Ai wi 

Wi \ 2 w 2 



J 



V 



W n _2 A n _l W„_i 



such that rifc=i w k 7^ 0. The Jacobi matrix can be decomposed into the product 

J = GJG (27) 

where G = diag(7i, 72, ... , 7 n ) is a diagonal matrix and J is a Jacobi matrix with all 
units on the neighboring parallels to the diagonal, 



Ai 1 

1 A 2 1 



J 



\ 



A„_i 1 

1 v 
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Hence det(J) = (IlLi7fc) det(J), and one can employ formula (I26p to evaluate det( J) 
(in the case of odd dimension). In more detail, one can put 



fc-i fe-i 

72fc-l = 



— , 72fe = wi I I , k = 1, 2, 3, ... . 



Alternatively, the sequence {7fc}fc =1 is defined recursively by 71 = 1, jk+i = Wk/lk- 
Furthermore, X\ = Xk/lk- With this choice, f l2"T|) is clearly true. 

Next we aim to derive a formula for the characteristic function of a Jacobi matrix 
with an antisymmetric diagonal. Suppose A = (A_<f, A_d_|_i, A_d+2, . . . , A d ) G C M+1 and 
A_fc = — Afc for — d < < d; in particular, Ao = 0. We consider the Jacobi matrix 
K = diag(A) + wE + + wE^. Let us denote, temporarily, by S the diagonal matrix with 
alternating signs on the diagonal, S = diag(l, —1, 1, . . . , 1), and by Q the permutation 
matrix with the entries Qj.fc = $j+k,a for — d < j,k < d. The commutation relations 

SQKQS = -K, S 2 = Q 2 = I, 

imply 

det(K-zI) = det(SQ(K - zI)QS) = - det(K + zl). 

Hence the characteristic function of K is an odd polynomial in the variable z. This 
can be also seen from the explicit formula ( 128]) derived below. 

Proposition 25. Suppose d £ N, w E C, A G C 2d+1 and A_fc = —Xk for k = 
—d, —d + 1, — d + 2, . . . , d. Then 

- - det ( diag(A) + wE + + wE_ - zl) (28) 



vfc=l 

d / d 



\x + z X d + Z / 

w W \ / W It) 



Aj+i + 2: A d + z 



Proof. This is a particular case of f l26|) where one has to set yk = Xk — z for k > 0, 
?/o = — = — (A_fc + z) for k < 0. To complete the proof it suffices to verify that 

s - — — u 



z(Ai — z) \X 2 — z X d — zJ \Xi + z X 2 +z 'Xd + z 
w 2 w w w \ u; u> 



5 1 — ,t — — d 



z (Ai + z) \ Ai - 2 ' A 2 — z ' " ' ' X d — z J V A 2 + z ' " ' ' A rf + z 

p[ Vfc=i k~ z ) vVi- 2 Ad-^y + 2 Ad + ^ 

To this end, one can apply (|S]), with n = d, Uk = w/(Xk — z), Vk = w/{Xk + z). Note 
that Uj — Vj = 2zUjVj/w. □ 
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Zero always belongs to spectrum of the Jacobi matrix K for the characteristic 
function is odd. Moreover, as is well known and as it simply follows from the analysis 
of the eigenvalue equation, if w 7^ then to every eigenvalue of K there belongs exactly 
one linearly independent eigenvector. 

Proposition 26. Suppose w G C, A G C 2d+1 , \- k = — X k for —d < k < d, and 
Tlfc=i 7^ 0. Then the vector v G C 2d+1 , v T = (9_ d , 9_ d+ i, 9_ d+2 , . . . , 9 d ), with the 
entries 



0* = (-l)VM If A, ]3 — . - for k = 0,1,2,..., d, (29) 




6-k = {— l) fc ^fc for —d < k < d, belongs to the kernel of the Jacobi matrix 
diag(A) + wE + + wE_. In particular, 9 = A X A 2 . . . \d$(w/ \i,w/ \ 2 , ■ ■ ■ , w/Xa), 
9 d = (—l) d w d , and sow^O. 

Remark. Clearly, formulas (121?]) can be extended to the case Ylk=i Afc = as well 
provided one makes the obvious cancellations. 

Proof. One has to show that 

w9 k -i + \ k 9 k + w9 k+ i = 0, k — — d + 1, — d+ 2, . . . , d — 1, 

and \- d 9- d + w9-a+i — 0, w9 d ^i + \ d 9 d = 0. Owing to the symmetries A_^ = —\ k , 
9- k = {—l) k 9 k , it suffices to verify the equalities only for indices < k < d. This can 
be readily carried out using the explicit formulas ( 129]) and the rule ([3]). □ 



5 Jacobi matrices with a linear diagonal 

Finally we focus on finite-dimensional Jacobi matrices of odd dimension whose diag- 
onal depends linearly on the index and whose parallels to the diagonal are constant. 
Without loss of generality one can assume that the diagonal equals 
{-d, -d+ 1, -d + 2, . . . , d), d G Z+. For w G C put 

K = diag(-d, -d + 1, -d + 2, . . . , d), K(w) = K + wE + + wE_. 

Concerning the characteristic function = det(K(w) — z), we know that this is an 
odd function. Put 

Xred{z) = det(K(w) - z). 

Hence Xre&{ z ) is an even polynomial of degree 2d. Further, denote by 
{e_d, e_d + i, e_d+2, • • • , e d } the standard basis in C 2d+1 . 

Suppose w^O. Let us consider a family of column vectors x s>n G C 2d+1 depending 
on the parameters s, n G Z and defined by 

x ^ n — {3( s + d, n),3(s + d- 1, n),3(s + d - 2, n), . . . - d, n)). 
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From the fact that the matrix 3 obeys (1181) . f)19p one derives that 

Ws,n G Z, K(w)x S:n = sx Stn — wZ(s + d + l,n)e_d — w$(s — d — l,n)ed- 

Put 

v s — ^s,s+d+l; S G Z. 

Recalling that Z(m, m) = 3(— m, m) = one has 

K{w)v s = s v s — w%(s — d — l,s + d + l)ed- (30) 
Remark 27. Putting s = one gets K(w)vo = 0, and so vo spans the kernel of K{w). 
Lemma 28. For every I = —d, —d + 1, — d + 2, . . . , d, one has 

W d+l (d + S )\ (£-s)\ Vse6i + s P an { e «+l' e ^+2; • • • ; e 4- 

In particular, 

d ( — -\\d+s 

„_,), '■• (3D 

s=—d 

Consequently, V = {v-d, V-d+i, • • • , ^d} a frasis m C 2a!+1 . 

Proof. One has to show that 



s=— <i 

Note that for any a G C and n G Z + , 



n\ f a + k 
n 



D-i)"C)( o : fc )=o.-o.u -i. b-d'CX 

fc=0 \ / \ / fc=0 \ / \ 

Using these equalities and fT24|) one can readily show, more generally, that 

/ T~i \7T7i rr + s, n + s) = for m, n G Z, m < n < m + d + 

z — ' (a + s)! — s)! 



•1)". 



and 



g d (d + X- 3)l ^ + s,m + d + e + s + l) = w- 



■d-e 



This proves the lemma. □ 
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Denote by K{w) the matrix of K(w) in the basis V introduced in Lemma 
Let a, b 6 C M+1 be the column vectors defined by a T = (a_d, a-d+i, a~d+2, ■ ■ ■ a^), 
& T = {P-d, P-d+i, P-d+2, ■ ■ ■ 0d), 

oc s = 3(s — d—1, s+d+1), (3 S = — -— -7, s = — d, — d+1, -d+2, . . . , d. (32) 

(a + sy {a — sy. 

Note that 

a- s = ~a s , (3-s = (3 S - (33) 
The former equality follows from ([22]) and (TT9]) . From fl30|) and fl3Tl) one deduces that 

K(w) = K - ba T . (34) 

Note, however, that the components of the vectors a and b depend on w, too, though 
not indicated in the notation. 

According to ( |34"|) . K{w) differs from the diagonal matrix K by a rank-one cor- 
rection. This form is suitable for various computations. Particularly, one can express 
the resolvent of K(w) explicitly, 

(K(w) - z)- 1 = (K - z)~ l + i _ ^ _ (K - zyViKo - z)-\ 

The equality holds for any z6C such that z ^ spec{K } = {— d, —d+1, —d+2, . . . , d} 
and 1 — a T (Ko — z)~ l b ^ 0. Clearly, this set of excluded values of z is finite. 

Let us proceed to derivation of a formula for the characteristic function of K(w). 
Proposition [231 is applicable to K{w) and so 

Below we derive a more convenient formula for Xred(^)- 
Lemma 29. One has 

d 

s\ {2d- 2s + 1)! 

and 



^ ((rf- g )!) 2 (2rf- g + l)! 2s 



fc=0 

/or n = 1, 2, . . . , d. 
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Proof. Let us first verify the formula for Xred(O). From (1351) it follows that 

x-W = (<«)'»(«,, f....,H) , + 2t^(f) , »( J ^ T 



2 

WW W ' 



By Proposition [T31 



Xred(O) = 7r 2 W M+2 Y d+1 (2w) 2 ^J {2w) 2 + 2 £ J 3 {2w) 2 J + 0( W M+2 logH) . 

Further we need some basic facts concerning Bessel functions; see, for instance, [TJ 
Chp. 9]. Recall that 

oo 

J (z) 2 + 2^J,(^) 2 = 1. 

Hence 

Xred (0) = ix 2 w 2d+2 Y d+l {2w) 2 + 0(w 2d+2 \og{w)) 

J2^l w A + 0(w 2d+2 \o g (w)). 



Note that Xrod(O) is a polynomial in the variable w of degree 2d, and so 

A -A (d- fc)!(d-s + fc)! 2s 

2.2. ^ «, . 

o— n I — n \ ' 



s=0 fc=0 

Using the identity 



E(d- k)\ (d — s + k)\ ,, , ,,x2 v-^ (d - k\ /. 

k=0 v 7 fe=0 v 7 v 

((d-s)!) 2 (2d-s + l)! 



d — k\ / d — s + k 
d — s J V d — s 



s! (2d — 2s + 1)! 

one arrives at fl36|) . 

To show ([37]) one can make use of f[M|) . One has 

_ (-\) d+l 
Xred(z) = det( K(w) -z) = i '- det(K - z) det(/ - (K - zYV) . 

z z v ' 

Note that det(7 + ba T ) = 1 + a T b. Hence, in view of (j3"3"D . 

xUz) = fl(k 2 - z 2 ) (i-£^l)= f[(k 2 - z 2 ) ( 1 - 2 £ ^ 



s — z I ^ ^ \ * — ' s 2 — £ 2 

fc=l \ s=—d / fc=l \ s=l 



18 



Using ( 132)) one gets 



_ V d 



Xred(n) = -2nf3 n a n \[ (k 2 - n 2 ) = ± '- w 2d+1 3(n - d - 1, n + d + 1). 

k=l 

k^n 

Formula (|37|) then follows from ( 124|) . □ 
Proposition 30. For every d G Z + one has 

xUz) = £ 7 + f]> 2 - z 2 ). (38) 

s=0 ^ ' k=l 

Proof. Since Xred(-z) is an even polynomial in z of degree 2d it is enough to check 
that the RHS of (13"81 coincides, for z = 0, 1, 2, ... , d, with Xred(0), Xred(l), Xred(2), . . . , 
Xred(^)- With the knowledge of values (13T)j) and (13T1) . this is a matter of straightforward 
computation. □ 

Remark 31. Using (I3"8"j) it is not difficult to check that formula (13"T|) is valid for any 
neN, including n > d (the summation index runs from 1 to min{n — l,d}). 

Remark 32. If w G R, w ^ 0, then the spectrum of the Jacobi matrix K(w) is real and 
simple, and formula (138 j) implies that the interval [—1,1] contains no other eigenvalue 
except of 0. 

Eigenvectors of K(w) can be expressed in terms of the function too. Suppose 
w^fl. Let us introduce the vector-valued function x(z) G C 2d+1 depending on z G C, 

u / \ j i, iYz + d+ 1) ^( w w w \ , 

With the aid of (EJ) one derives the equality 

_ 2d T(z + d + 1) ^( ww w 



(K {W ) - #) = J^, - 7TT , . . . , — d j e, (39) 

Remark 33. According to ( TT2]) . 

s=0 v ' j=s 

Hence £k( z ) is a polynomial in z of degree d + k. In particular, = 1, and so 

x{z) ^ 0. 

Proposition 34. One has 



ww w w , ( 

z — d'z — d + l'z — d + 2'' ' z + d 



If w G C, u> 7^ 0, i/ien /or every eigenvalue A G spec(iT(ty)) ; x(A) zs an eigenvector 
corresponding to A. 
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Proof. Denote by P(z) the RHS of (HP]) . By ([22]), if P(A) = then x(X) is an eigenvec- 
tor of K(w). Thus it suffices to verify fHUj) . The both sides depend on w polynomially 
and so it is enough to prove the equality for w G R \ {0}. Note that P(z) is a poly- 
nomial in z of degree 2c? + 1, and the coefficient standing at z 2d+l equals —1. The set 
of roots of P(z) is contained in spec(-fr(u>)). One can show that P(z) has no multiple 
roots. In fact, suppose P(X) = P'(X) = for some A G R. From f[3"9l one deduces that 
(K(w) — A)x(A) = 0, (K(w) — A)x'(A) = x(A) (here x'{z) is the derivative of x(z)). 
Hence 

(K(w) - A) 2 x(A) = (K(w) - A)V(A) = 0. 

Note that x'(A) ^ 0, and x'(X) differs from a multiple of x(A) for C,-d{ z ) — 1- This 
contradicts the fact, however, that the spectrum of K(w) is simple. One concludes 
that the set of roots of P(z) coincides with spec(i^(w)). Necessarily, P(z) is equal to 
the characteristic function of K(w). □ 

Remark 35. With the aid of (HOI) one can rederive equality f!3T|) . For 1 < n < d, a 
straightforward computation gives 

x ( n ) = (-lf+n w 2d+l^ d _ n + 1 ^ d + n+1 y 
Equality (J3TD then follows from fl24j) . 
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